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Abstract 

. A schematic model of over-damped motion is presented which permits one 

\ to calculate the mean first passage time for nuclear fission. Its asymptotic value 

CD ' may exceed considerably the lifetime suggested by Kramers rate formula, which 

applies only to very special, favorable potentials and temperatures. The addi- 
^r~i \ tional time obtained in the more general case is seen to allow for a considerable 

' increment in the emission of light particles. 



PACS numbers: 24.75. +i, 24.60.Dr, 24.10.Pa, 24.60.-k 



Q . 

^ ■ Fission experiments are commonly analyzed with the help of statistical codes 

which are based on simple rate formulas for the processes of fission and emission of 
light particles. In the early 80'ies an excess of neutrons was observed over that for 
^ \ which the fission rate is simply estimated by the Bohr Wheeler formula Ff = Fbw 

' (see e.g. the review articles PQE] with references to original work). An improvement 

was seen in replacing Fbw by the Fk of Kramers jS] in which the fission rate formula 
gets reduced by friction, the more the larger the dissipation strength. Additional 
possibilities for enhancing the relative emission probability Fn/Ff of light particles 
like neutrons were attributed to two effects which seem to arise in a time dependent 
description: (i) Starting the dynamics of fission at some time zero, it takes a finite time 
for the current across the barrier to reach thestationary value from which Kramers 
derived his formula, (ii) To this stationary current a finite time lapse Tssc may be 
associated for the motion from the saddle point down to scission. Often feature (i) 
is interpreted as a delay of fission during which particles may be emitted on top of 
the number given by Fh/Fk. Likewise, it is believed that also the neutrons emitted 
during Tssc are not accounted for by this Fh/Fk. 

A review of these features and of their practical applications can be found in jSj. 
It can be said that interesting consequences have been deduced in this way, both for 
the value of the dissipation strength as well as for its variation with temperature, 
see e.g. jUIH]- More recently, however, the question has been raised as to whether 



1 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 


0.4 0.8 1.2 1.6 2.0 

Q 

Figure 1: The mean first passage time t{Q = Qcx) = Tmfpt(Qa Qcx) of (jT} (nor- 
malized to tlie tk of 0) for a cubic potential and different temperatures, defined in 
units of E'b- 




Kramers' original rate formula itself accounts for realistic situations in fission [H]. For 
under-damped motion, modification becomes necessary whenever the inertia changes 
from the minimum to the barrier. Moreover, it has been argued that any temperature 
dependence of the pre-factor must not only be attributed to friction, but also to the 
inertia and in particular to the stiffnesses of the potential at its extrema. In the 
interpretation of fission decay as a sequence of three subsequent steps (minimum- 
saddle, motion across saddle, saddle-scission) has been re-examined with the help of 
the concept of the "mean first passage time" (MFPT). Restricting to over-damped 
motion such an analysis can be performed in analytic fashion, simply because for this 
case an analytic formula for the time r^fpt exists. It reads 
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and is valid if any coordinate dependence of friction 7 and temperature T are dis- 
carded, details may be found in [S]. Here, Qa is meant to represent that minimum 
of the potential V{Q) which is associated to the "ground state deformation" and Qex 
stands for the "exit point". In this sense the r^fptiQa. Qcx) determines the average 
time the system spends in the interval from Qa to Qex- It is calculated for a situation 
where the system, after starting at Qa sharp, does not return to this interval once it 
has crossed the point Qcx, which is then referred to as an "absorbing barrier" . Typical 
for fission, for Q —>■ — cxd the V{Q) is assumed to rise to plus infinity, and hence acts 
as a "reflecting barrier". 

As will be demonstrated again below, the r^fptiQa. Qex) tends to a constant 
as soon as the exit point is sufficiently far to the right of the potential barrier. This 
constant, which henceforth shall simply be called Tmfpt, becomes identical to the 
inverse of Kramers' rate r^fpt = tr = ^/^k, whenever the usual conditions are fulfilled 
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which underly Kramers' derivation. Recalhng that we are deahng with over-damped 
motion, this tk is given by 
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exp(Eb/T) , 



(2) 



Tk = 



where and Cb are the stiffnesses at the potential minimum and barrier, respectively. 
In jS] this fact is proven by applying the saddle point approximation to formula (jT)). 
For this it is important to have exactly two saddle points, those corresponding to one 
minimum and one barrier. Notice that the saddle point approximation requires one 
to replace the barrier by an inverted oscillator, which indeed was also assumed to 
hold true by Kramers in his famous work. 

As a typical case we show in Fig|T]the results of calculations of rjnfpt(<5a <5cx)/tk 
for the same cubic potential as used in [7]. It may be specified by its first derivative 
to be given by the form 



with the barrier height -Eb = V{Qh) — V^(<5a) to be 8 MeV with the extrema to lie at 
— 0.2 and — 0.8. It is observed that the asymptotic ratio rmfpt/rx becomes 

close to unity, indeed, if only the parameter temperature over barrier height becomes 

small enough. However, even for this case of exactly two well pronounced extrema, 

deviations from unity are clearly visible at larger temperatures. 

This situation becomes more dramatic as soon as the potential shows additional 

structure. This will now be demonstrated using a schematic potential of fifth order 

in Q. Again, V{Q) will be fixed by its first derivative. 



with i?b5 Qa and Qb unchanged. The remaining two parameters Qc and Qd niay be 
used to specify structure of the potential beyond the barrier. In FigsO they have 
been chosen to be identical to one another, Qc = Qd = Qs, with their values fixed 
such that the height of the then existing shoulder takes on the values specified in 
the figure captions. 

In all cases the calculation of r^iptiQa Qcx) was performed up to regions of the 
exit point Qex where the stationary value is reached. It is seen that this asymptotic 
regime is not very far away from the one where the potential is assumed to have 
additional structure. This is so even for the example shown in the lower right corner 
of Fig|21 There a potential is taken with a shoulder in a region which lies 20 MeV 
below the first minimum, or — 2.5£^b in terms of the barrier height. For heavy nuclei 
this may thus be said to correspond to the scission region after which the fragments 
separate. The ratio rnifpt(Qa Qcx)/tk shown in the figures is calculated for the tk of 
(j21). Evidently friction drops out but the stiffnesses Ca and Cb from © remain. They 
are taken to be those of the individual potentials for which the Tmfpt is computed. 
These results exhibit clearly the mistake one makes if only Kramers' rate formula 
is used to estimate the time the system stays together. Rather, the considerable 



v'iQ) (X (g-ga)(g-Qb) , 



(3) 



V'iQ) oc (g - ga) (g - gb) (g - a) (g - Qd) , 



(4) 



3 





2- 



V(Q)/E 



0.5- 



X(Q=Qe.)/XK 




0.2- 
1.0- 



0.5 



1.0 q1.5 



2.0 2.5 



V(Q)/Ej . 


^(Q=Qex)/^K 








-T/E.= 0.1 


- \ 0.2 


^ 0.5— 


- \ 1.0 







0.5 1.0 1.5 r\ 2.0 2.5 3.0 



Figure 2: Same as in Fig^but for potentials having shoulders at some Qs of heights 
Vs/ Ey, relative to the barrier: V^/-Eb = 0.75 top left, 0.5 top right, bottom left, —2.5 
bottom right. 



overshoot of Tmfpt over tk indicates that much more time is available for light particles 
to be emitted before scission. Suppose we look at neutrons. Whenever, their average 
width Fn may be used to calculate their multiplicity per fission event from Fn/Ff the 
enhancement of this number over that given by Fn/Fx is determined by the ratio 

Tmfpt /t"K, viz 

Fn Tjnfpt , , 

T^ = T^ • 

if i-K Tk 

For the potentials chosen here this enhancement may become quite large. 

To get some feeling for absolute values of this extra available time we estimated 
the pre-factor 7/ a/ Ca|Cb| of Q following suggestions given in [H]. We simply re- 
placed the geometric mean A/Ca | Cb | by the C which there in sect.in.B.5 appears in 
the relaxation time Tcow for over-damped collective motion. Rather than use the for- 
mulas given in that section we simply take the results shown in Fig. 4 (by the dashed 
line). Between temperatures of one and four MeV this Tcou shows an almost linear 
dependence in T such that one may write 

7 3 5^ / n \ 

^ - -7 + 7^ Ut77 • (6) 



Putting this estimate into the formula given in (0) for tk{T) one obtains the curve 
shown in FigOl The strong temperature dependence reflects the exponential function 
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Figure 3: The MFPT corresponding to Kramers' estimate as function of temperature 
in units of the barrier height, see text. 



exp(£'b/T). As the estimate © seizes to be vahd above T ~ 4 MeV the curve 
should not be taken too seriously above TjE\^ ~ 0.5. For such a T the tk is about 
200 /i/MeV large. As the rmfpt/rx typically is about 1.5 the additional time increment 
= Tmfpt — Tk takes on the sizable value of roughly 100 h/MeV, and, hence, is at 
least as large ClS db typical transient time. 

Finally, in Fig0 we look at the case of the potential having a second minimum 
and maximum. As was to be expected the effect is even larger than before. We 
would like to remark, however, that this example should be taken with some caution. 
Commonly, such a double humped barrier comes about because of shell effects. In 
the range of temperatures considered here the latter may be considered to be quite 
weak if not already washed out completely. After all our study is concerned with 
over-damped motion. According to j^] (see also [H]) nuclear collective motion may be 
expected to become over-damped only above temperatures of about T = 2 MeV. 

Our results may be summarized as follows. One of the main issues has been to 
corroborate features suggested before in jTj, and to some extent already in 0. In 
essence they imply the following two issues: (i) For situations for which transport 
equations like those of Kramers or Smoluchowski (or the corresponding Langevin 
equations) may be applied to analyze fission experiments, it may be inadequate to 
simply use Kramers' famous rate formula. Deviations from that may originate in 
various reasons, for instance in transport coefficients varying with shape, see [Bj. Here, 
we concentrated on properties of the potential restricting ourselves to over-damped 
motion and the model case of constant friction, (ii) For this model we have been able 
to demonstrate that there is considerable room for increasing the time the fissioning 
system stays together without having to rely on concepts which are meaningful only 
within a time dependent picture. Whereas results obtained within the latter may 
depend crucially on initial conditions this is not the case for the MFPT This Tmfpt 
represents the average time it takes for the system to start at the potential minimum 
and to make its motion all the way out to scission. It includes relaxation processes 
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Figure 4: Like in FigI21 but for a potential having an additional minimum of depth 
Vc = — 0.25£'b and a second maximum of height = 0.25 -Eb- 



around the first minimum as well as the sliding down from saddle to scission. The way 
it is derived |H] implies a proper incorporation of averages over the statistics which are 
to be associated with a process underlying fluctuating forces. Whereas for Kramers' 
model case the Tmfpt is nothing else but his inverse rate, this is no longer true for larger 
temperatures and, in particular, not for potentials of more complicated structure. We 
have been able to demonstrate that the latter may lead to a considerable prolongation 
of the time the systems spends before it scissions, allowing in this way for emission 
of light particles on top of those given by T^/T^. Of course, further work will be 
necessary to clarify the relevance of this feature with respect to real situations. For 
such studies not only more realistic potentials have to be used, one should also try 
to generalize formula (0) to include a coordinate dependent friction coefficient. The 
ultimate goal should be to be able to use such a formula in a statistical code where 
one may account for temperatures changing during the process. 
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